To appear in ApJ letters October 10, 2003 



Volatile-Rich Earth-Mass Planets in the Habitable Zone 

Marc J. Kuchner 1 

ABSTRACT 

A small planet is not necessarily a terrestrial planet. Planets that form beyond 
the snow line with too little mass to seed rapid gas accretion (< 10 M@) should 
be rich in volatile ices like H 2 and NH 3 . Some of these planets should migrate 
inward by interacting with a circumstellar disk or with other planets. Such 
objects can retain their volatiles for billions of years or longer at ~ 1 AU as their 
thick steam atmospheres undergo slow hydrodynamic escape. These objects could 
appear in future surveys for extrasolar Earth analogs. 

Subject headings: astrobiology — planetary systems 



1. INTRODUCTION 

Several upcoming space missions 2 aim to detect extrasolar planets with masses of order 
1 M0. The Space Interferometry Mission (SIM) should be able to detect extrasolar planets 
down to ~ 3 M0 astrometrically by observing the reflex orbital motions of their host stars 
(Ford & Tremaine 2003). The Kepler and COROT missions are designed to photometrically 
detect planets with the same radius as the Earth by observing the dimming of their host stars 
during planetary transits. Eventually, missions like the Terrestrial Planet Finder (TPF) and 
Darwin should directly detect extrasolar planets around nearby stars and take their spectra, 
searching for planets as faint as the Earth in either reflected light (Kuchner & Spergel 2003) 
or thermal emission (Woolf 2003). 

Earth-mass planets offer particular biological interest if they orbit in the habitable zone, 
the range of circumstellar radii where stable liquid water can exist on the surface of a planet 
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that is in thermal equilibrium with stellar radiation (Huang 1959; Kasting et al. 1993). In the 
solar system, this zone covers roughly the region between the orbits of Venus (0.72 AU) and 
Mars (1.52 AU). Natural places to look for life-bearing extrasolar planets are the habitable 
zones of nearby stars. 

Discussions of extrasolar planets often quietly assume that any object with mass ~ 
orbiting in a star's habitable zone will be terrestrial, i.e, composed mostly of silicates and 
iron-peak elements — like the Earth. However, we suggest that the habitable zones of nearby 
stars could harbor other similar-looking beasts. Planets composed substantially of volatiles 
can conceivably form with masses up to ~ IOMq and migrate inward to the habitable zone. 
We estimate that Earth-sized volatile-rich planets can retain their volatiles for billions of 
years at ~ 1 AU, protected by thick atmospheres that slowly undergo hydrodynamic escape. 

2. FORMATION AND MIGRATION 

The snow line in a protoplanetary disk (Hayashi 1981) occurs where the temperature 
of the nebular gas drops below the sublimation temperature for water (170 K in a vacuum). 
In popular planet-formation scenarios (e.g. Pollak et al. 1996; Wuchterl 2000), planets form 
much faster beyond the snow line than in the warmer center of the disk because outside the 
snow line, frozen volatiles provide more solid building material. Protoplanets beyond the 
snow line more massive than a threshold mass roughly equal to the presumed core masses 
of Jupiter and Saturn rapidly accrete gas, becoming gas giants (but see also Boss 2002). 
Protoplanets beyond the snow line less massive than this threshold grow slowly, and when 
the gaseous disk has dispersed they remain small balls of ice and rock, like Pluto and the 
icy moons of the outer planets. 

Icy planets and moons that form near the snow line should contain mostly H2O and 
silicates. Those that form farther out can acquire NH 3 and CH 4 as hydrates, together 
with CO, N 2 , and CO2. Beyond the snow line, the solar nebula does not have time to reach 
chemical equilibrium, making it hard to predict the relative ratios of these compounds (Prinn 
1993). Icy bodies more massive than ~ 0.02 form with enough accretional energy that 
they begin life with their ice mostly melted (Lunine & Stevenson 1982; Stevenson et al. 
1986). These objects must freeze again in the outer solar system, though the details of this 
process remain uncertain. 

Jupiter's moon Ganymede has a mass of only 0.02 but it makes a good reference 
point (Sohl et al. 2002). Ganymede has a diameter larger than that of Mercury (0.06 M®); 
its mean density suggests that it is roughly 60% silicate rock and metal and 40% volatiles. 
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Water ice appears to cover most of Ganymede's surface; this ice may top an ocean main- 
tained by radiogenic heating (Spohn & Schubert 2003). Objects even richer in water could 
possibly form near the snow line, where water may collect due to diffusion in the solar nebula 
(Stevenson & Lunine 1988). 

Interaction with a circumstellar disk can cause protoplanets to migrate inward from 
the snow line (Goldreich & Tremaine 1980; Lin et al. 1996). Ice planets could form with 
masses up to the threshold for rapid gas accretion (~ IOMq); low mass objects in this 
category undergo "type I" migration in a massive protoplanetary disk, which appears to 
occur on timescales of 10 5 (M/Mq)~ 1 years, where M is the mass of the planet (Ward 
1997). If the volatile-rich planet migrates after the ~ 10 7 year period of terrestrial planet 
formation at ~ 1 AU, it has a better chance of surviving near the star. Even after the outer 
protoplanetary disk has dissipated, interaction with other protoplanets can scatter an icy 
body into an orbit with small pericenter (Rasio & Ford 1996; Adams & Laughlin 2003), 
where further encounters, stellar tides (Mayor & Queloz 1995) or disk tides (Ward 1988) can 
circularize its orbit. Scattered Earth-mass planetesimals in the outer solar system may have 
generated the obliquity of Uranus and Neptune (Stern 1991) and the high inclination Kuiper 
Belt objects (Petit et al. 1999). 

A volatile planet may not even need to migrate to get to the habitable zone. Recent 
models of flared disks around young stars suggest that the snow line could be as close-in as 
1 AU in a disk around a solar- mass star (Sasselov & Lecar 2000). Asteroids orbiting between 
2.6 and 3.5 AU (most C-type asteroids) contain minerals that appear to have undergone 
aqueous alteration, suggesting that they contained water ice when they formed (see the 
review by Barucci et al. 1995). A snowline located closer-in than 2.6 AU would naturally 
explain the existence of icy primitive asteroids. 



A water planet arriving at ~ 1 AU should quickly develop a steam atmosphere of 100 
bars or more, like the atmosphere of the accreting Earth (Abe & Matsui 1985). Such an 
atmosphere can easily become a runaway greenhouse when accretional impacts deliver an 
energy flux > the solar constant to its base. But in the absence of such an energy source, the 
two-stream model of Matsui & Abe (1986) suggests that a steam atmosphere of > 30 bars 
simply provides stable greenhouse heating up to a temperature of 
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VAPORIZATION: A SELF-LIMITING PROCESS 
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at its base, where A is the ratio of the opacity at thermal wavelengths to the opacity at short 
wavelengths, u is the albedo for single-scattering at short wavelengths, and Tbb is the local 
blackbody temperature; T BB = 278 K (a/1 AU)- 1/2 (L*/L ) 1/4 , where a is the radius of the 
planet's orbit, and is the luminosity of the star. Since A is likely to be in the range of 
10-100 (A w 31 for Venus), and u is likely to be almost 1, T G « 2AT BB - 

When T G < 647 K, the critical temperature of water, water simply condenses out of 
a steam atmosphere that is in equilibrium with stellar radiation, leaving a hot ocean on 
the planet's surface (Matsui & Abe 1986; Zahnle et al. 1988). In this situation, the planet's 
vaporization halts, and the planet can survive indefinitely as long as it retains its atmosphere. 
This condition holds where a > 1.03 AU(L*/L ) 1/2 , for A = 31, uo = 1. 

Closer in, the planet will slowly vaporize (or sublimate, if ice remains on its surface) at 
a rate set by the rate of energy transport through the base of the optically thick atmosphere, 
where the temperature will fall from Tq to 647 K at the surface of the boiling ocean. The 
atmosphere would be stable against convection since it is heated from above; downward 
energy transport in this stable layer will occur by slow diffusion of radiated thermal energy. 
The flux of energy transported by radiative diffusion through the lowest scale height of 
the atmosphere would be F di ff ps — (16/3)oT 3 AT/t r (T), where a is the Stefan-Boltzmann 
constant and tr(T) is the Rosseland mean optical depth of the layer. The temperature 
difference across the layer, AT will be AT = Tq — 647 K. The mean temperature, T, will be 
roughly (T G + 647 K)/2. 

The more massive the atmosphere, the slower the vaporization; this process is probably 
self-limiting. The vaporization rate is ~ A^r^Fdiff / L, where tq is the radius of the planet, and 
L 2.26 x 10 10 erg g _1 is the latent heat of vaporization of water. The optical depth of the 
atmosphere will increase as more volatiles sublimate, but when ps 10 5 (a/l Af/) _2 (L^/L Q ), 
the vaporization rate reaches ~ 10 -10 Mq yr^ 1 , i.e., it effectively halts. 

Suspended aerosols in the atmosphere could easily provide a continuum mass opacity of 
k — 0.1 cm 2 g _1 , like clouds in brown dwarfs (Cooper et al. 2003). The planet should have 
a radius similar to that of the zero temperature sphere described by Zapolsky & Salpeter 
(1969); a zero temperature sphere with mean atomic mass per atom of 13 AMU (half H 2 0, 
half Si0 2 ) has a radius of r ~ 1.4 R.0 (M/Mq) 1 / 3 . With these assumptions, the limiting 
optical depth amounts to a mass of gas 

M max = 2 x 10" 3 M (k/0.1 cm 2 g- 1 )- 1 (a/l AUy 2 {U/ L & ){M / M^f'\ (2) 

Though this thick atmosphere can protect the planet's volatiles from vaporization, the 
planet is vulnerable during the first billion years of the solar system's life, when the star 
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shines brightly in the EUV, and the planet suffers heavy bombardment by planetesimals. 
We will consider each of these processes in turn. 



4. EUV-DRIVEN ESCAPE 

Stellar EUV photons can raise the kinetic energy of gas molecules above the escape 
velocity, eroding the planet. For example, Pluto has probably lost a few km of surface ice 
during its lifetime to an EUV driven wind (see the review by Trafton et al. 1996). We can 
estimate the escape rate by estimating the radius, r±, where the atmosphere first becomes 
optically thick to EUV. This parameter sets the effective cross-section of the planet for 
absorbing stellar EUV and also the thickness of the layer which must conduct the deposited 
energy to the surface. We will neglect the vertical thickness of the opaque region of the 
atmosphere just above the surface; its scale height should be h s 59 km (M/M^) -2 / 3 
(assuming steam at 647 K), much less than the typical scale height of the upper layers. This 
approximation works as long as most of the planet remains liquid (or ice) so the atmosphere 
has only a small fraction of the planet's mass. 

Watson et al. (1981) and Hunten & Watson (1982) modeled the escape of Pluto's at- 
mosphere and the atmospheres of the primordial Earth and Venus using a one-dimensional 
hydrodynamic model based on the model for the solar wind by Parker (1963). This model 
estimates the effective value of r 1 for the case of maximum escape flux by assuming all the 
EUV energy is deposited in a thin layer at r\. This maximum escape flux provides a lower 
limit for the lifetime of the bound gas. 

Watson et al. (1981) define the variables A = GMm/kT r , (3 = S(GMm/kT$K ) and 
C = F esc (k 2 T /KoGMm), where m is the mass of a gas molecule. The base of the exosphere, 
which we will consider to be located at ro, has temperature T , which we take to be 2ATbb- 
The flux of escaping particles per steradian per second is F esc , and the effective EUV heating 
rate is S, in erg cm~ 2 s _1 , averaged over a sphere. The thermal conductivity of the gas is 
assumed to take the form k = k (T/T ) s . The maximum escaping flux satisfies (Watson et 
al. 1981): 

\\ = —. ^ r , (4) 

C (a - (p/0 1/2 ) 

where p — 2/(s + 1). We solved these equations numerically for ( as a function of f3 and Ao 
and integrated the mass loss rate, M = 47imF esc , over time, t, since the birth of the planetary 
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system. Since vaporization immediately replaces mass lost to the EUV-driven wind, M is 
the mass evolution of the planet. 

This hydrodynamic model applies when A > 10, i.e., the atmosphere is tightly bound. 
As A decreases, EUV penetrates deeper into the atmosphere, and r\ approaches r (see 
figure 3 of Watson et al. 1981). We assume that when A < 7, ri=r , meaning that the 
EUV penetrates all the way to where the atmosphere becomes opaque near the surface of 
the planet. In this case, F esc ks Sr^/(GMm). 

The thermosphere may have complicated photochemistry; EUV can dissociate any of the 
key atmospheric molecules (H 2 0, NH 3 , CH 4 , CO, N 2 , and C0 2 ). But we can conservatively 
estimate the loss rate of these compounds as the escape rate for monatomic H (m = 1.67 x 
10~ 24 g). Banks & Kockarts (1973) calculate that the thermal conductivity for H is k = 
379T o.69 erg cm -i s -i K -i between T = 200 and 2000 K. 

The effective value of the present day EUV flux at the Earth is S ~ 1 erg cm" 2 s" 1 
(Watson et al. 1981). Younger solar- mass stars can easily have 100 times as much EUV 
output as the sun. We assume that S = 1 erg cm -2 s" 1 (a/1 AU) at t = 4.5 billion years and 
thereafter. We assume that prior to that, the EUV output evolved as (t/4.5 x 10 9 years)~ ' 8 , 
up to S = 100 erg cm" 2 s" 1 (a/1 AU)" 2 following Zahnle & Walker (1982). 

Figure 1 shows the EUV survival time — the time for an all- volatile planet to lose half of 
its mass to an EUV-driven wind — as a function of initial mass for planets at a = 1, 0.3, 0.1, 
and 0.04 AU. The figure suggests that an Earth-mass volatile-rich planet as close as 0.3 AU 
from a sunlike star can easily retain its volatiles for the age of the solar system. Our model of 
EUV-driven escape assumes that most of that planet, say 90%, remains ice or liquid, so that 
the base of the thermosphere sits near the planet's surface. We estimate that this situation 
occurs only where a > 0.14 AU(M/M0)" 1 / 6 (L*/L Q ) 1 / 2 . The dashed sections of the curves 
in Figure 1 indicate the regime where this condition is not satisfied. 
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Fig. 1. — Survival time for EUV-driven atmospheric escape as a function of mass for volatile 
planets at 1 AU, 0.3 AU, 0.1 AU, and 0.04 AU. The curves become dashed where the 
atmosphere becomes > 10% of the planet's mass. The dotted line shows the age of the solar 
system. 
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Since volatile-rich Earth-mass planets slowly lose mass to a EUV-driven wind, these 
planets should be more common around younger stars and stars with lower time-averaged 
EUV luminosity. Nearby F-, G-, and K-type main-sequence stars have similar levels of EUV 
output, but nearby M-type main-sequence stars have systematically lower EUV luminosities 
than more massive nearby main-sequence stars (Wood et al. 1994). This difference might 
make M-type stars particularly hospitable to close-in, volatile-rich Earth-mass planets. 

5. BOMBARDMENT 

We assume that the volatile-rich planets migrate after the usual period of terrestrial 
planet accretion in habitable zone. But "heavy" bombardment by leftover planetesimals 
continued on the terrestrial planets in our solar system for roughly the first billion years of 
the solar system's life. Such a bombardment can erode the atmosphere of a planet (see the 
review by Ahrens 1993). 

A low-mass volatile-rich planet will have escape velocity v esc ~ 9 km s _1 (M/Mq) 1 / 3 . 
Ejecting a mass m e of material from the planet requires energy \m e v 2 esc . But the impactors 
available in the early solar system would rarely have impact velocities more than a few times 
v esc , so an impactor probably could not eject more than roughly its own mass from the 
planet's atmosphere, given the ~7% efficient ejection process found in numerical models of 
large impacts on the Earth (Newman et al. 1999). Moreover, only impactors with diameter 
greater than the scale height of the atmosphere impart a substantial fraction of their energy 
to the atmosphere (O'Keefe & Ahrens 1982); smaller impactors just add mass to the planet. 
A scale height of 59 km corresponds to a projectile with mass ~4 x 10 17 kg. 

We appropriate a description provided by Melosh & Vickery (1989) for the flux of 
impactors on Mars and the Moon based on crater densities. The number of projectiles per 
year of mass greater than m p that fall on a planet of effective radius r e ff is 

N cum = 0.25 yr- 1 (r e// /%) 2 (l + 2300e-'/ 2 - 2xl ° 8 K/lkg)" ' 47 , (5) 

where we have neglected gravitational focusing — a factor of a few when scaling from Mars 
to a 1 planet. If we consider only impactors with masses between 10 15 kg and 1% of the 
mass of the moon on a planet with r e // = Rq, Equation 5 implies a total mass of impactors 
of less than 0.003 This total mass of impactors poses no hazard whatsoever, even to 

the survival of an 0.1 volatile-rich planet. 

However, the above impactor distribution fails at the high mass end; it does not account 
for the unusual impactors (Tremaine 1993) that have masses comparable to a terrestrial 
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planet. Such giant impactors probably created the Earth's moon (e.g., Ward & Cameron 
1978; Canup & Asphaug 2001) and may have spun up Mars (Dones & Tremaine 1993) and 
robbed Mercury of its low-density material (Benz et al. 1988). Conceivably, a giant impactor 
with high enough velocity could remove all the volatiles from a 1 volatile planet. 

6. CONCLUSIONS 

The proximity of the snow line to 1 AU in accretion disk models suggests that volatile- 
rich planets can easily migrate into the habitable zones of solar-type stars. As long as they 
migrate after terrestrial planet formation and manage to avoid giant impacts, we estimate 
that volatile-rich planets as small as 0.1 Mq can retain most of their volatiles for the lifetime 
of the solar system as close as 0.3 AU from a solar-type star. Closer-in objects may also retain 
volatiles, but our analysis fails at high heating rates, where the atmosphere's mass becomes 
substantial. Exterior to a > 1.03 AU (L i(i /L ) 1 / 2 , a water planet can exist in equilibrium 
with stellar radiation, protected by a thick steam atmosphere. Inside this radius, the planet 
slowly vaporizes, but we estimate that only a mass of water given by Equation 2 can vaporize 
in 10 10 years. Conceivably, both the Earth and Mars could once have been volatile-rich, but 
both have probably suffered large impacts. Volatile-rich planets should be more common in 
the habitable zones of young stars and M-type stars. 

Volatile-rich Earth-mass planets could appear in future surveys for extrasolar Earth 
analogs. These objects should maintain deep layers of liquid water, melted by accretional 
energy and radionucleides. At a given mass, a volatile-rich planet should have a slightly 
larger radius and therefore a greater luminosity than a terrestrial planet. Finding such 
objects would test core-accretion scenarios for giant planet formation. A small planet is not 
necessarily a terrestrial planet. 

Thanks to Ann Bragg, Ulyana Dyudina, Mike Lecar, Dimitar Sasselov, Sarah Stewart, 
and Wes Traub for discussions. Thanks also to Alain Leger for sharing with me some of his 
unpublished work on ocean planets. This work was performed in part under contract with 
the Jet Propulsion Laboratory (JPL) through the Michelson Fellowship program funded by 
NASA as an element of the Planet Finder Program. JPL is managed for NASA by the 
California Institute of Technology. 
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